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Abstract- This  paper  presents  a  generalized  approach  to  derive  the  dyadic  Green’s  func¬ 
tion  of  arbitrary  multilayer  structures  in  planar,  cylindrical  and  spherical  coordinate 
systems.  It  is  based  on  a  full-wave  equivalent-circuit  representation  that  makes  it  possible 
to  apply  simple  network  analysis  techniques.  The  planar  stratified  structures  may  be  lat¬ 
erally  bounded  or  open,  the  cylinders  or  spheres  circumferentially  closed  or  limited  sec¬ 
tors.  The  procedure  can  be  extended  to  other  coordinate  systems  and  anisotropic  materi¬ 
als  by  means  of  a  generalized  transmission-line  approach.  For  a  demonstration  of  the 
applicability  of  this  method,  the  derivation  of  the  dyadic  Green’s  functions  in  spectral 
domain  of  some  commonly  used  cylindrical  and  spherical  microstrip  structures  is  pre¬ 
sented. 


I.  INTRODUCTION 

Apart  from  several  other  applications,  the  dyadic  Green’s  function  is  often  found  as 
kernel  in  integral-equation  technique  in  combination  with  the  method  of  moments 
(MoM)  to  solve  the  boundary  value  problem  of  microstrip  antennas.  These  may  con¬ 
sist  of  several  layers  dielectric  with  metallizations  in  the  interfaces  and  can  have  pla¬ 
nar,  cylindrical,  spherical  or  any  other  shape.  A  microstrip  structure  in  a  general  coor¬ 
dinate  system  is  depicted  in  Fig.  1. 


Fig.  1:  Conformal  multilayer  microstrip  stmcture  in  curvilinear  coordinate  system. 

The  derivation  of  the  dyadic  Green’s  function  can  be  done  by  the  use  of  a  full-wave 
equivalent  circuit  and  was  already  derived  for  certain  coordinate  systems  [1],  [2],  [3]. 
The  general  description  in  curvilinear  coordinates  has  been  presented  in  [4].  The  pre¬ 
sent  paper  briefly  describes  the  general  concept  behind  and  shows  two  applications  of 
microstrip  antennas  in  cylindrical  and  spherical  coordinates. 
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II.  ANALYSIS 


Starting  from  the  wave  equations  in  different  coordinate  systems  the  relation  between 
the  tangential  field  components  in  spectral  domain  on  both  sides  of  a  layer  can  be  de¬ 
rived  [l]-[4] 
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The  derivation  can  be  done  either  by  following  a  general  curvilinear  description  [4]  or 
in  addition  to  this,  the  derivation  of  above  equation  was  done  via  the  vector  potentials 
for  the  spherical  case  in  [3].  Here  the  tangential  field  components  have  been  derived 
from  the  vector  potentials  and  can  be  described  by  the  coefficients  A,  5,  C,  D  and  the 
according  Hankel  functions  using  four  equations.  The  coefficients  are  eliminated  by 
setting  up  the  equations  on  top  and  bottom  of  a  layer  k  from  which  eq.  (1)  is  obtained. 
Furthermore  the  boundary  conditions  at  the  interfaces  are  taken  into  account  by 

(2) 

On  bottom  and  top  of  the  whole  structure  a  relation  between  electric  and  magnetic  field 
components  is  set  up 


H-0,/2  _  E0,w  5 
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where  the  admittance  Y0  n  is  determined  by  either  having  open  or  closed  boundaries 
below  or  above  the  multilayered  structure. 

Eqs.  (l)-(3)  are  the  starting  point  for  the  general  description  of  multilayered  structures 
by  the  use  of  a  full-wave  equivalent  circuit.  Fig.  2  shows  the  equivalent  circuit,  contain¬ 
ing  the  above  equations. 
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Fig.  2:  Full-wave  equivalent  representation  of  a  conformal  multilayer  stmcture,  where 
is  the  first  of  the  curvilinear  coordinates. 


Finally  the  dyadic  Green’s  function  in  spectral  domain  is  obtained  [1] 

E  =  G  J  (4) 

Now  the  solution  via  integral  equation  technique  in  combination  with  spectral  domain 
approach  can  be  used  to  solve  the  boundary  value  problem  of  microstrip  antennas.  Ex- 
emplarily  this  is  presented  for  microstrip  antennas  in  two  coordinate  systems,  the  cy¬ 
lindrical  and  the  spherical  one. 


III.  APPLICATION 

Fig.  3  shows  the  first  application.  A  single-layer  microstrip  antenna  structure  is  de¬ 
picted  in  cylindrical  coordinates.  Furthermore  the  related  equivalent  circuit  is  shown. 
The  dyadic  Green’s  function  for  the  single-layer  microstrip  structures  is  obtained  as 
[5] 


G  =  (Zr1F1+F2)-1  .  (5) 

The  matrix  components  can  be  obtained  as  described  in  [5]  or  the  general  formalism  in 
[4]  can  be  used.  The  procedure  is  the  same,  only  in  [5]  the  cylindrical  coordinates  are 
used  right  from  the  beginning,  whereas  in  [4]  the  derivation  is  done  for  general  curvi¬ 
linear  coordinates  and  later  on  these  are  replaced  by  cylindrical  ones. 


a)  b) 

Fig.  3:  a)  Single-layer  microstrip  antenna  structure  in  cylindrical  coordinates  and  b)  re¬ 
lated  equivalent  circuit. 


From  the  solution  of  the  integral  equation  the  unknown  currents  are  obtained  and  pa¬ 
rameters  like  input  impedance  (Fig.  4)  or  far  field  patterns  can  be  derived  [5]. 


Fig.  4:  Input  impedance  of  a  direct  fed  ^-oriented  microstrip  antenna  [5].  fo=2. 2-2  A  GHz 
(clockwise),  p0= 30  mm,  d\=p\-p0=OJ%l  mm,  sY i=2.33.  Simulation  (o),  measurement 

(n). 


The  general  procedure  is  now  applied  to  the  case  of  a  microstrip  patch  antenna  on  a 
circumferentially  closed  spherical  structure  of  one  dielectric  layer  covering  a  metallic 
sphere  as  shown  in  Fig.  5.  The  corresponding  equivalent  circuit  is  the  same  as  shown  in 
Fig.  3b). 


1  Metallizations 

2.  Substrates 


Fig.  5:  Microstrip  antennas  on  a  circumferentially  closed  spherical  structure. 


The  Green’s  function  for  this  structure  is  expressed  as  already  given  in  Eq.  (5).  The 
matrix  elements  are  now  defined  by  the  spherical  coordinates,  where  the  sub¬ 
matrices  Y2  ,  Zx  and  Vx  of  the  hybrid  matrix  can  be  found  in  [3]. 

The  calculated  input  impedance  of  the  cp-polarized  patch  is  compared  with  commercial 
software  as  shown  in  Fig.  6. 


Fig. 6:  Input  impedance  of  phi-polarized  patch  ri=60  mm,  r0=58.42  mm,  £ri=2.52  quasi- 
rectangular  cp-polarized  patch  with  sizes  of  w9=25  mm,  we=40  mm. 

IV.  CONCLUSIONS 

A  general  procedure  for  the  derivation  of  the  dyadic  Green’s  function  of  arbitrary  mul¬ 
tilayer  structures  has  been  presented.  Two  applications  of  single-layer  microstrip  an¬ 
tennas  in  cylindrical  and  spherical  coordinates  have  shown  the  applicability  of  the  pre¬ 
sented  approach. 
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